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^^ ! Tlie antiferromagnetic Heisenberg Hamiltonian is investigated on a truncated 

tetraliedron, which is a ciosed 12 site system. We find that the ground state has many 
similarities to that of Cqo • We study 2- and 4-spin correlations in the classical ground 
state of the truncated tetrahedron and calculate the same correlations in the exact 
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Q> , S=l/2 ground state. We find that the classical correlations survive for a range of 



2 ! bond strengths in the Heisenberg Hamiltonian and that one can construct a good trial 

__i , wavefunction based on the classical ground state. This suggests that the correlations 

Q , present in the classical ground state of Cqq also survive in the exact ground state of 

O 

that system, for a range of bond strengths about the physically relevant J2 ~ Ji. 

rS i We calculate the momentum-space correlation function S{q), which is measurable 

by neutron scattering, for both C12 and Cqq. We also calculate correlations at finite 

temperature. 
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I. INTRODUCTION 

We have found that an antiferromagnetic Heisenberg Hamiltonian on Cqq has an exotic 
ground state with nontrivial topology. ||l| The Heisenberg Hamiltonian arises from a tight- 
binding Hamiltonian describing the half-filled vr orbitals at the carbon sites, assuming that 
there is strong Coulomb repulsion between electrons in these orbitals. The repulsion may 
be either long or short range. The exotic classical {S = cxd) ground state has nontrivial 
topology, with a skyrmion number 7. (The vector field n that is normal to the spin plane 
of the pentagons covers the sphere 7 times.) Cgo is, however, a spin 1/2 system, and 
quantum corrections are expected to be important. Estimates using perturbation theory 
indicate that the corrections are about 60% to the ground state energy and about 40% to 
the magnetization per site. Given the size of these corrections, there is concern that the 
perturbative results may be misleading. 

The truncated tetrahedron C12 is system that is similar to Cqq, but is small enough that 
the S = 1/2 problem can be solved exactly. The classical ground state for C12 also has 
nontrivial topology, with a skyrmion number 1. Here we investigate 2- and 4-spin correlations 
in the exact ground state of the S = 1/2 Heisenberg antiferromagnet on Cu, and compare 
to correlations in the classical ground state. We find that the classical correlations survive 
in the exact ground state for a wide range of the parameter J2/J1, the ratio of the two bond 
strengths in the system. We also construct an S = 1/2 variational wavefunction based on a 
coherent state representation of the classical ground state spin configuration, and show that 
it has a large overlap with the exact ground state. 

These results suggest that the correlations found in the classical ground state of Ceo also 
survive for S = 1/2 for a range of bond strengths. There is expected to be no phase transition 
in the exact ground state as a function of on-site repulsion U in the Hubbard model on the 
Ceo lattice. [|I],^ The spin correlations for the large U limit Heisenberg antiferromagnet 
are therefore expected to be present, although somewhat reduced in magnitude, for the 
physically relevant case of intermediate U. 



II. CLASSICAL APPROXIMATION 

As in the case of Cqo, all carbon atoms on C12 are equivalent, but there are two inequiva- 
lent types of bonds. The magnetic exchange constant for two neighboring sites on the same 
triangle is Ji = t\/U^ where ti is the hopping matrix element and U is the Hubbard inter- 
action. The constant connecting a site on one triangle with a nearest neighbor on another 
triangle is J 2 = t\/U. The Hamiltonian is 

t nt 

H = Ji Yl ^j-^k + J2 Yl ^3 ■ '^fc' (1) 

<hk> <j,k> 

where the first sum is over the 12 triangle bonds, and the second over the 6 non-triangle 
bonds. We take Ji =1 and investigate the properties of this Heisenberg Hamiltonian on C12 
as a function of J2- We first treat the Hamiltonian classically, so that a in Eq. (|l|) is a 
classical unit vector. This is the same as the S —>■ 00 limit. We discuss the properties of the 
exact quantum {S = 1/2) ground state in the next section. 

The truncated tetrahedron has 4 triangular faces and 4 hexagonal faces. An antifer- 
romagnet is frustrated on the triangular faces. For an isolated triangle, the ground state 
energy is Ji cos(27r/3) = — Ji/2 per bond, which is smaller than the — Ji per bond that 
an unfrustrated system would have. If the 12 spins on the truncated tetrahedron are able 
to find an arrangement so that the different triangles do not further frustrate each other, 
the ground state energy will reach the lower bound Eh = 12 Ji cos(27r/3) — 6 J2. The lower 
bound is an energy — Ji/2 per triangle bond, and —J2 per non-triangle bond. The classical 
ground state configuration was found numerically by minimizing the energy over the spin 
variables {6'j,0j}, i = 1, 12 (see Fig. 1). We find that C12 does attain the lower bound Ej,, 
but only in a topologically nontrivial spin arrangement. In comparison. Ceo also attains its 
lower bound, but other fullerenes C20, Cjo, and Cs4 do not. ||l[| 

The ground state configuration for C12 is the same for all positive Ji and J2- In addition 
to the global rotational symmetry of the ground state, there is a discrete parity symmetry 
whereby each spin at -^ —cji. The state related to a ground state Gi by parity cannot be 



reached by a global rotation on G*i. Although the ground state spin configuration is non 
co-planar, it may be easily visualized. There is a global rotation such that the spins on 
each triangle lie in the plane of that triangle. The skyrmion number s for a set of vectors Vi 
defined at points i on the surface of a sphere is given by 

where a is the solid angle subtended by three vectors and the sum is over a set of fundamental 
triangles that cover the sphere. The skyrmion number for C12 and Ceo is calculated by 
introducing "twist" bond variables Ui on the triangle bonds for the vectors Vi, with 

hi ~ aj X ak, (3) 

and one finds that both of the parity-related ground states of Cu have skyrmion number 
1. (The skyrmion number cannot be defined directly on the spin variables because the solid 
angle in Eq. (|^) is undefined for many sets of 3 neighboring spins.) The hi defined by any 
pair of spins on the same triangle are equal. The spin configuration may be chosen so that 
the hi point along the normals to the triangles as shown in Fig. (1). The sum of all hi is 
zero. The net magnetic moment in the ground state is also zero. 

In the presence of an applied magnetic field, Cu has a second order metamagnetic tran- 
sition at which dM/dh is discontinuous. The skyrmion number is 1 both above and below 
the transition. In contrast, C20 and Cgo, both of whose ground states have skyrmion number 
7, have a first order metamagnetic transition at which M is discontinuous. 

III. SPIN 1/2 

For 5* = 1/2, the a matrices in Eq. (|1|) represent 2x2 Pauli spin matrices. (Our 
convention is to use a matrices, which are a factor of 2 larger than s matrices.) The 
Heisenberg Hamiltonian on the truncated tetrahedron was solved in each Sz sector separately. 
H The dimension of the S^ = Q sector is 924 and that of the S, = 1 sector is 792, so that 



the total spin 5* = subspace has dimension 924 — 792 = 132. The lowest energy eigenvalue 
is present in the 5*^ = sector and absent from the Sz = 1 sector, so the ground state 
is a singlet. We calculate various correlation functions as a function of J2/Ji- In Cqq, J2 
is expected to be nearly equal to, but slightly larger than Ji, because the corresponding 
non-pentagon bond is about 3% shorter than the pentagon bond. [^ C12 has not been 
synthesized, but we consider J2 ~ Ji to be the region of physical interest. 

First we consider the quantum two-spin correlations Qij =< a-i ■ (fj >. There are five 
inequivalent gi/s on the truncated tetrahedron, which are represented by gu, gu, gi5, gis 
and gig, where the subscripts refer to the site labeling in Fig. (1). The gij^s, with the 
exception of gig, are plotted in Fig. (2). (The classical gig is zero, and the quantum 
Igigl < 0.1 for all J2.) The functions gu and gi4 are the triangle and non-triangle bond 
energies divided by the bond strengths. One sees that as J2 -^ 0, gi2 approaches -1, the 
value on an isolated triangle. Surprisingly, as J2 ^ 0, the correlation function gi4 remains 
large and negative. When J2 = 0, there is no correlation between spins on different triangles. 
Then < ai ■ 0^4 > is averaged over a manifold of degenerate ground states and is equal to 
zero. For J2 small, however, the expectation is over the unique ground state and is nonzero. 
In the opposite limit of large J2, gu approaches -3 and all other gij fall to zero. This reflects 
the fact that in this limit, the wavefunction becomes a product of singlet valence bonds on 
the 6 non-triangle bonds, denoted \VB >. The values of gi^ and gis, the third and fourth 
nearest neighbor correlation functions, both peak at Ji ~ J2. 

We now calculate correlations between the spin plane normals fii in the exact ground 
state. This will involve four-spin correlation functions. As we pointed out in the discussion 
of the classical ground state, the sum of the spin plane normals Ylt=i ^i ^^ ^^e four triangles 
is zero, so that classically 

hi ■ {n2 + hs + Hi) = -1. (4) 

These correlations may be measured by 

G =< hi ■ (n2 + h3 + h^) >, (5) 



where 

na = i^i X ^j + Oj X Ok + (Tfc X aj)/3. (6) 

The spins on triangle a are numbered counterclockwise {i,j,k). (The Ua in Eq. (Q) are 
normalized to unit length, but the tIq, in Eqs. (^§]) are not.) 

The function G measures the relative orientation of the spin-plane normals, and involves 
correlations between spins over the whole system. It is plotted in Fig. (3). As J2 increases 
from zero, \G\ increases slowly to its maximum at J2/J1 ~ 0.3. For large J2, \G\ approaches 
zero as the wavefunction approaches \VB > and all correlations other than those for the 
singlet on nearest neighbor non-triangle bonds vanish. The four-point function G falls off 
at somewhat smaller values of J2 than the two-point functions in Fig. (2). The quantum 
correlation functions will be compared with classical and other correlation functions below. 
Note that the skyrmion number itself is not a good quantum number for the quantum 
system. 

Singh and Narayanan ^ have emphasised that in identifying a nonvanishing expectation 
value of an operator in the ground state as evidence for magnetic order, it is important to 
show that the expectation value is considerably larger at zero than at infinite temperature. 
Note in this regard that G and all of the Qij are traceless operators and thus have zero 
expectation value at T = 00. 

A classical ground state has each spin pointing in a fixed direction. The most straight- 
forward translation of a classical ground state into a spin 1/2 wavefunction is to have each 
of the quantum spins point in the same fixed directions. The wavefunction for the first spin 
is the two-component spinor |ci > that is the ground state oi Hi = —hi-ai, where hi points 
in the classical spin direction. These are coherent states. The full wavefunction is 

\C >= \Ci > |C2 > . . . |Ci2 > . (7) 

Since the classical ground state configuration is independent of J2/J1, so is the trial wave- 
function \C > (and other trial wavefunctions discussed below). The expectation of the spin 
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correlation functions in the trial wavefunction \C > and in the exact groundstate {ipex > for 
Ji = J2 are shown in Table 1. The exact ground state has large 2-spin and 4-spin correla- 
tions, even at large separations. The correlation functions Qij {i 7^ j) and G are identical in 
\C > and in the classical ground state. The state \C > does reasonably well with all of the 
correlation functions, including G. The largest error is the underestimate of gi4, which has 
some singlet character in \ipex >, but can have only Ising character in \G >. 

Any global rotation of a classical ground state gives another ground state. Thus a better 
trial state would be to sum \G > over all global rotations, or equivalently to project it onto 
the singlet subspace. As mentioned in Section II, the parity operation a, — > — a^ results 
in another set of ground states that cannot be reached by rotations. We therefore define 
states l^ps > and I?/'! > as the (normalized) singlet projections of the parity related coherent 
states. An additional related state is found by taking a linear combination of these, 

|^5>=ai|^^>+a2|^|>, (8) 

where ai and 02 are complex numbers with |aip + |a2p = 1. We determine Oi and 02 in Eq. 
(^) by maximizing | < ipsl'^ex > P- (Virtually identical results are obtained by minimizing 
the energy.) The amplitudes ai and 02 are found to have equal magnitudes and a phase 
difference that is almost independent of J2/ Ji- In Fig- (4) we plot the norm of the overlap 
squared of these states and of \VB > with the exact ground state {ipex >• One sees that 
I < i^sli'ex > P and I < i'sl'ipex > P are both peaked near Ji = J2, and that l'?/'^ > is a good 
approximation to the ground state when J2 ~ Ji. \VB > becomes exact for J2 '^ Ji, but is 
already failing in the region of physical interest, J2 ~ Ji. 

The correlation functions of the trial states \VB >, {ipg >, and {ips > are given in Table 
1. (1^1 > has the same correlation functions as 1^^ >, and is not shown.) The state \VB > 
gives no indication of either the magnitude or sign of any spin correlation function other than 
gi4. As the variational wavefunction is improved in the sequence \G >^ \ips >^ {ips >, 
the correlation functions generally approach those of the exact ground state. Convergence 
is best at short distances and worst for the long distance gis- By examining Figs. (2-3) and 



Table 1, one sees that the trial state {ips > gives a reasonable approximation to all of the 
correlation functions not just for J2 ~ Ji, but also as J2 -^ 0. The agreement is remarkably 
good considering that zero-point spin-wave fluctuations have not been added to any of the 
trial states. 

We now compare the energy of trial states and that obtained from perturbation theory 
to the exact ground state energy. High energy contributions from non-singlet states are 
removed by projecting \C > onto the singlet subspace. The energy of {ips > is thus lower 
than that of |C >. The energy of \i/js > is -8.803 Ji -12.325 J2, which is an upper bound 
on the exact ground state energy. For Ji = J2 = 1, this energy is -21.128, which is to 
be compared with an exact ground state energy of -22.804. From perturbation theory, the 
leading correction to the ground state energy is 

ae=i:'^aw^_ (9) 
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eo-ej 



where |0 >= \C > and the sum is over intermediate states \j > with two adjacent spin flips. 

(The ground state is determined by stability against single spin flips.) For Cu this is 

72 27 72 

AE=-Q^---^ -. (10) 

Ji 2 (2J2 + Ji) ^ ' 

The perturbed ground state energy is E^"^^ = Ed + AE, with Ed = — 6(Ji + J2). This 
expression clearly overestimates the correction to the classical ground state energy when 
J2 > Ji, but it does fairly well for J2 ^ Ji- For J^ = J2 = 1, E^^) ^ -22.500, which should 
be compared with the exact ground state energy of -22.804. 

For Cqo, we expect the trial state \VB > to be less successful than for C12 in competing 
with the other trial states. The reason is that Cu contains highly frustrated triangle bonds, 
which can be sacriflced at relatively low cost by the formation of singlets on the non-triangle 
bonds. A pentagon is less frustrated than a triangle, with a classical energy per bond of -.809, 
as compared to a triangle's -0.5 and an unfrustrated bond's -1. It is therefore reaonable to 
expect that the correlations found in the classical ground state of Cqo, which with skyrmion 
number 7 are more exotic than those in the truncated tetrahedron, should survive for a 
larger range of J2 values. 



Neutron scattering experiments can measure S{q), the Fourier transform of the spin-spin 
correlation function. We calculate S{q), which is S averaged over all relative orientations of 
g and the molecule. This would be appropriate for a sample containing molecules of random 
orientations. Neglecting form factors, 

'S'(g) =J2<^i'^J> sm{qrij)/qrij, (11) 

where Vij is the cartesian distance between sites i and j. We assume that the triangle and 
non-triangle bonds lengths are equal. The calculated S(q) for lipex >, li^s >, \VB >, and 
the coherent state \C > are shown in Fig. (5). (S(q) for |C > is not zero at q=0 because 
|C > is not a singlet.) If the triangle and non-triangle bond lengths differ by very little, as 
is the case for the pentagon and non-pentagon bonds in Ceo, it is difficult to separate the 
contributions to S(q) from gi2 and gu. The longer distance spin correlations, present in all 
of the wavefunctions except \VB >, cause higher frequency oscillations in S{q). This leads, 
for example, to a shoulder at qr 12/11 ~ 2.7, which is strongest in |C > but still visible in the 
exact ground state. 

For Cqo, we do not know [ipex > and show S{q) for the coherent state \C > obtained 
from the classical ground state, and for the state \VB > with singlets on the non-pentagon 
bonds. There is a great deal more structure in S(q) for \C > than for \VB >. These features 
will, however, be smoothed out somewhat for the exact ground state wavefunction for Cqq 
as they were for Cu- 

Finite temperature: We have calculated the C12 correlation functions Qij at finite tem- 
perature, using the exact eigenfunctions with Boltzmann weight. As expected, all Qij for 
distinct spins approach zero as T -^ 00. Two of the correlation functions, however, are 
nonmonotonic. The correlation gi2 between spins on the same triangle increases in magni- 
tude by 6% with increasing T before decreasing. The correlation gig, which is quite small at 
zero temperature, actually changes sign and increases in magnitude by a factor of 16 before 
decreasing to zero. 

The finite temperature properties can help to interpret recent Quantum Monte Carlo 
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studies on the Ceo Hubbard model with 60 electrons by Scalettar et al. [^ They find that for 
U /t = 8, the spin correlations on Cqq resemble those of the classical ground state at short 
and intermediate distances, but decay to near zero at the largest separations. A crucial 
question is whether this long distance decay is a property of the quantum ground state, or 
whether it is caused by the finite temperature of the simulation. For Cu, we calculated 
the exact S = 1/2 Heisenberg partition function for J corresponding to U/t = 8 and the 
lowest temperature Ti = t/5 reached. The temperature Ti = 0.582A, where A is the gap 
to the first excited state. The first excited state is, however, 9-fold degenerate counting 
spatial and spin degrees of freedom, and higher excited states are fairly close. At T = Ti, 
the system has a probability 0.19 to be in the quantum ground state, and the correlation 
function ^fig is 48% of its zero temperature value. Ceo should require a lower temperature 
than Ci2 for the long distance spin correlations to reach their zero temperature values. // 
the Ceo ground state had spin correlations extending across the molecule, one would expect 
the Ceo excitation energies, which correspond to magnons at allowed nonzero wavevector. 



to be smaller than those of Cu by roughly a factor of two (~ W 12/60). By this estimate. 
Ceo at a given temperature would have similar properties to C12 at twice that temperature. 
At T = 2Ti, C12 has a probability 0.028 to be in the ground state, and a long distance 
correlation function that is reduced by a factor of 6 from its zero temperature value. This 
suggests that the decay of spin correlations at large distances observed in the Quantum 
Monte Carlo studies may be a finite temperature effect. 

IV. SUMMARY AND CONCLUSION 

The classical Heisenberg antiferromagnet on the truncated tetrahedron has an exotic 
ground state with nontrivial topology and nonzero spin correlations up to the longest dis- 
tances in the system. To see whether quantum fiuctuations for S = 1/2 destroy these corre- 
lations, we exactly diagonalized the S* = 1/2 problem. We find that the classical ground state 
gives a very good indication of the exact quantum 2-spin and 4-spin correlation functions 
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for Ji K. J2. The classical ground state is also used to generate trial S = 1/2 wavefunctions 
with good overlap against the ground state for Ji k, Jj- Although the wavefunction \VB > 
with singlet valence bonds on the non-triangle bonds is the correct one for Ji ^ J2, this 
state gives no inkling of the long range correlations that develop for the Ji ~ J2 region of 
physical interest. 

This work has implications for the magnetic properties of Ceo, in which an exotic classical 
ground state has been found. The results for C12 suggest that the S = 1/2 quantum 
fluctuations will not destroy the classical spin correlations in Ceo. The quantum fluctuations 
may play a smaller role in Ceo, because the spin arrangement is less frustrated than that of 
Ci2- There is also concern that spin correlations in Ceo can be destroyed because U/t ~ 4 
is not large enough. We have argued that the exact quantum ground state of the Hubbard 
model in Ceo is unlikely to have a phase transition as a function of U/t, and thus that the 
large U spin correlations, somewhat reduced in magnitude, should be observable in real Ceo. 
|T|] Evidence for these magnetic correlations may be found in inelastic neutron scattering 
experiments. It may also be possible to infer the spin correlations by doing NMR on a Ceo 
molecule with two ^^C nuclei replaced by ^^C. It remains to investigate how doping affects 
the spin correlations, and to what extent spin correlations affect the properties of doped Ceo 
crystals, such as KsCqq. |§| 
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TABLES 
TABLE L Correlation functions for the trial wavefunctions and for the exact ground state 

{ipex > with J2 = Jl- 





\VB> 


\c> 


l#> 


\i^s> 


IV'e^ > 


912 





-0.5 


-0.7745 


-0.7336 


-0.7295 


9lA 


-3 


-LO 


-L5890 


-2.0541 


-2.3416 


915 





0.5 


0.4761 


0.5131 


0.4453 


918 





-0.5 


-0.7745 


-0.7323 


-0.4916 


919 








-0.1088 


-0.0418 


0.0013 


G 





-0.75 


-0.4879 


-0.5548 


-0.5947 
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FIGURES 
FIG. 1. Spin configuration on a truncated tetrahedron. A spin pointing up out of the plane has 

a larger head, and one pointing into the plane has a larger tail. For the global rotation shown, all 

spins lie in the physical plane of their triangle, which is the plane of the paper for triangle (123). 

The truncated tetrahedron has been flattened by stretching the back hexagon. The spin directions 

are shown unchanged. 

FIG. 2. Two-spin correlation functions gij =< ai ■ (Tj >, 5'i2(dot-dash line), gi^ldash line), 
gi5(full line), and (7i8(dot line), as a function of J2/ Ji- 

FIG. 3. Four-spin correlation function G, defined in Eq. (||), as a function of J2I J\- 

FIG. 4. Overlap squared of the trial wavefunctions with the exact ground state, | < V'S'lV'ea; > P 
(full line), I < V'sl^ex > P = I < V'llV'ex > P (dot-dash line), and | < VB\ti)ex > P (dot line). 

FIG. 5. The angle-averaged spin correlation function S{q) for C12 is plotted as a function of 
qru/n for \ipex >(full line), 1-05 >(dot line), \VB >(dot-dash line), and \C >(dash line). 

FIG. 6. The angle-averaged spin correlation function S{q) for Cqo is plotted as a function of 
qri2/iT for \C >(full line) and \VB >(dash line). 



14 



